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Multiple Timescale Dynamics in Economic Production Networks 
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We describe a new complex system model of an evolving production economy. This model is 
the simplest we can envisage which incorporates the new observation that the rate of an economic 
production process depends only on the minimum of its supplies of inputs. We describe how this 
condition gives rise to a new type of complex multiple timescale dynamical evolution through a 
novel type of bifurcation we call a trapping bifurcation, which is also shown to be one cause of non- 
equilibrium economic behaviour. Such dynamics is an example of meta-level coupling which may 
also arise in other fields such as cellular organization as a network of molecular machines. 

PACS numbers: 89.65.-s, 89.65.Gh, 89.75.Fb, 05.45.-a, 82.39.Rt, 05.45.Xt 



Introduction. Recently physicists have been devot- 
ing a great deal of attention to economic and financial 
phenomena |l]|. Economic dynamics is easily observed to 
be far from equilibrium where periodic recessions, un- 
employment and unstable prices occur persistently. An 
understanding of the origins of this behaviour from the 
viewpoint of complex dynamical systems theory would be 
very valuable. Our model @], which we hope takes a step 
in this direction, is based on von-Neumann 's3] neoclassi- 
cal model of economic production and catalytic chemical 
reaction network dynamics 

The original von Neumann model (VNM) of economic 
production assumes that each good is produced jointly 
with certain others, in an analogous way to a chemical re- 
action. A production process is the operation which con- 
verts one bundle of goods, including capital equipment, 
into another bundle of goods, including the capital equip- 
ment. Capital goods therefore function approximately 
like catalysts in chemical reactions, reformed at the end of 
the reaction in amounts conserved in the reaction. Con- 
sumption of goods is assumed to take place only through 
production processes, including life necessities consumed 
by workers, and all income is reinvested in production. 
Therefore the VNM is defined by a fixed input matrix 
and a fixed output matrix, representing, the fixed sto- 
chiometric ratios of input products and output products 
for each process. An example of such an economic net- 
work is shown in Fig^ Each process is assumed to have 
unit time duration, longer processes being broken down 
into several processes with intermediate products. The 
VNM is defined as a static equilibrium model describing 
relationships between the variables which must hold at 
equilibrium. Equilibrium is a state of 'balanced growth' 
where prices are constant. There are no dynamics defined 
by the model which might describe out of equilibrium or 
approach to equilibrium behaviour. 

However it is rarely the case that economic processes 
are in equilibrium. For example consider a simple bak- 
ery with input: (dough, baker, oven) and output: (bread, 
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FIG. 1: An Economy. Red lines are production flows, blue 
lines flows to and from the environment (including other coun- 
tries) and green lines are catalytic effects. Some processes are 
autocatalytic. 



baker, oven), so that the baker and the oven are catalysts. 
We observe that the rate of bread production, and eco- 
nomic processes in general, depends on the quantity of 
the minimum of its input supplies and not on the quan- 
tities of its other supplies. I.e. the baker working at full 
pace can only fill ovens at a certain rate. Similarly em- 
ploying more bakers will not increase bread production 
if the ovens are already full. This is in direct contrast 
with a chemical reaction obeying the law of mass action, 
where an increase in concentration of any of its input 
species will increase the rate of reaction. 

This minimum condition leads directly to multiple 
timescales]J$ . Consider a factory assembling bicycles 
from wheels and frames at a given rate. Suppose for some 
reason the rate of supply of wheels increases so that their 
price drops. To maximise efficiency the factory will try to 
use all its funds to expand production, buying the extra 
wheels and demanding more frames. A factory with large 
funds may demand frames so strongly that the the frame- 
makers respond to it and increase their supply rate. In 
this way the wheel-makers and frame-makers production 
rates may become synchronised. A factory with small 
funds on the other hand may only weakly influence the 
frame supply rate. Its bicycle production rate will be 
determined by the supply rate of whichever is the most 
expensive, and may switch between the rates as price lev- 
els change. The factory will have inefficient periods with 
surpluses of frames or wheels, but to avoid this by not 
using all its funds would still not be maximally efficient. 



Here we show how such multiple timescales arise 'en- 
dogenously' in a new model of an economy of coupled 
processes which is different from the VNM in that we ex- 
plicitly take into account the way a process' production 
rate depends on the minimum of its input supplies. 

Model. Although our model ignores many details, we 
believe it captures the essential characteristics of a evolv- 
ing production-marketing network, in the simplest way 
possible. The system is defined by a fixed number of 
processes Kj and products Kj with input stochiometric 
ratios ay and ouput stochiometric ratios fry, where i la- 
bels the processes and j the products. Each process has 
supplies of its products Sij(t), it also has a dynamical 
fitness called funds Fi(t), measured in dollars, which rep- 
resents the 'size' or 'intensity' of the process. The model 
can be conveniently broken down into four simple parts. 

(i) Processing. Here all processes simultaneously man- 
ufacture goods so that their supplies after processing, 
denoted with '*', are given by, 

Stf (*) = Sij(t) + {kj - aiJ )Min k (S lk {t)/a tk ) (1) 

where the processing rate is given by Mink{Sik{t) / a^) 
where Min k (x k ) denotes the minimum over x k . 

(ii) Product revaluation. The value, or price pj (t) , of a 
product is recalculated after processing and depends on 
the overall supply available in the whole economy and 
on how much all the processes want it. To measure how 
much the processes want the product we consider how 
much they will pay for it and each process therefore di- 
vides its funds Fi(t) into demands Dij(t) for its input 
products. The new product value Pj(t) is then given by 



Pj(t) = 



Total Demand JJi D%j(t) + Dj 



Total Supply 



(2) 



where Sj and Dj are possible external supplies and de- 
mands coming from, say, another country. 

There are many ways a process may allocate its funds 
into demands. If we consider a feedback from the price 
the funds should be divided according to aijPj(t — 1) 
since this is process i's best estimate of the per unit cost 
of input j depending on the most recent known price 
Pj(t — 1). Alternatively the process may simply divide its 
funds according to the ratios . Here we will consider, 

aijPjit - 1) 



D^t) = Fi(t)<n 



(1 - 7K,) (3) 



where the parameter < 7 < 1 quantifies the feedback 
from the price. In fact we will show that the interest- 
ing behaviour is found independently of 7, the important 
points being -Dy(t) oc Fi(t) Vj and £ • Dij(t) = Fi(t). 

(iii) Product resharing. Processes obtain new supplies 
for processing simply according to 



S ij (t+l)=D ij (t)/ Pj (t) 



(4) 



so that processes with relatively large funds obtain larger 
proportions of the available supply. 
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FIG. 2: Three examples of multi-timescale time series pro- 
duced from 15 process, 10 product networks, (7 = 0.9). 



(iv) Process revaluation. The new process funds are, 

F(t + i) = J2s; J (t)p 3 (t) (5) 

j 

so that processes with larger manufactured supplies ob- 
tain greater proportions of the available funds, and in 
this way the system becomes an evolving network. 

Results. These equations have the novel property that 
they switch their form depending on which product is the 
minimum at any time. Synergetic effects then appear 
because the dynamics itself selects the equations which 
determine its own future evolution. We will show that 
this strong non-linearity leads to the new type of very 
complex multi-timescale dynamics illustrated in FigEl 
focusing on 7 = since the analysis is much simpler. 

We first show that even the simplest single process, 
in a fixed environment, will, under general conditions, 
oscillate periodically due to a novel switching mechanism 
and furthermore this has interesting economic relevance. 

This Kj = 1, Kj = 3 system has one input product I, 
one catalyst C, and one output product O, reacting with 
non-zero stochiometric ratios ai = ac = be = bo = 1/2 
with 7 = 0. It has 3 variables, its input supplies Si(t), 
catalyst supplies Sc{t) and funds F(t), where we drop 
the process label i. There are two phases: (X) excess 
input, Si(t) > Sc(t) and (Y) excess catalyst, Si(t) < 
S c {t) 7 such that EqU S* 7 (t + 1) = l/2F(t)/ Pl (t) and 
S c (t+1) = l/2F(t)/p c (t) becomes, 



(X) Sj(t+1) = 
(Y) 



l/2F(t) + D c j 

s l{ t + i)= l ' 2F ^ 



(X,Y) S c (t + l) 



l/2F(i) + Dj 
l/2F(t)(S c (t)+S* c ) 
l/2F(t) + D c 



(6) 
(7) 
(8) 
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where we have omitted the (X) and (Y) F(t) equations 
for simplicity and where Sf, Sq, Sq and Df, Dq, Dq are 
fixed external supplies and demands of the three products 
respectively. Here, since 7 = and ai — ac, (X) and (Y) 
switch whenever the prices pi(t) and pc{t) cross. 

To understand the novel dynamical features of this sys- 
tem it is sufficient to consider the interaction of the phase 
(X) and phase (Y) fixed points. Denoting these two fixed 
points by superscript ' x ' v ' we find that, although the fixed 
point positions are different, the fixed point prices in both 
phases are the same and given by, 



Pi 



x.y 

Po 



Sj ~\~ Sq 



x,y _ ~C 

Pc — ~nr- 



(9) 



Therefore a process in equilibrium acts to make its in- 
put and output product prices equal 6] . This interesting 
non-trivial relationship has the observable consequence 
that for example scaling both D Q and Sq by the same 
constant factor will change the price of both O and I, 
although one would naively expect that such a scaling 
should leave the prices unchanged. 

However as shown in FigOJa), which shows a typical 
time series variation in the Sj(i), Sc(t) plane, the pro- 
cess will not in general be at a fixed point, but oscillates 
periodically. Although neither phase (X) nor (Y) has a 
periodic attractor on its own, a limit cycle is created by 
a novel trapped switching mechanism which arises due 
to the minimum dynamics as explained in Fig|2a). As 
shown the ratio of the fixed point values p = (Pq V /pp v ) 
acts as a trapping bifurcation parameter since when p > f 
the switching state is obtained but when p < 1 the sys- 
tem remains at the node without oscillating. 

The switching state is plausible economically. In the 
switching state we find the fund fixed points are such 
that F x — pF y > F y . Therefore when the system is in 
the excess catalyst phase (Y) its funds F(t) decrease and 
it acts to reduce Sc(t) (EqJHJ. However the decrease in 
Sc(t), since it is not the minimum supply, has no effect 
on the processing rate, and Si(t) (EqQl is therefore not 
affected by the change in Sc(t). In the excess supply 
phase (X) on the other hand, the funds F(t) increase and 
the process acts to increase Sc(t) (EqJSJ which increases 
the processing rate and decreases Si (t) (EqUjJ producing 
the oscillation. Furthermore as explained in Fig[3fa)(c) 
the switching frequency decreases as p increases. This is 
to be expected economically. The price fixed points are 
determined by the external environment. When p^, y » 
p X j V the process demands will have much less affect on 
Pc{t) than pi(t). Therefore in the excess supply phase 
(X) Sc(t), and therefore the processing rate, will only 
increase slowly and Si(t) will be reduced only slowly, 
producing large amplitude slow oscillations. 

We therefore expect that even a single process in a 
fixed environment will have periodic dynamics, with pe- 
riods of (X) full employment (Si(t) > Sc(t)) interrupted 
by periods of (Y) unemployment, {Si(t) < Sc(t)). 

Through simulations we find Eqs|5]also seem to hold 
in the 7 > case. The same switching state appears, see 
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FIG. 3: (a) The switching state; creation of periodic or- 
bit from node and focus. The Si = Sc line (see key) di- 
vides the (X) and (Y) phases, whose fixed points are shown 
labelled. The fixed points are such that they lie on the 
Si — Sc(Pc V /pT V ) ' me i ( see key). Therefore only one can 
be compatible with its phase. Here, in the switching state 
shown, where Pc' v /pT V ) > 1 the (X) fixed point is compat- 
ible since Sf > Sq, but the (Y) fixed point is incompatible 
since Sj > Sq. When the trajectory is in the (Y) phase 
it is attracted to the stable node (S c ,Sj), which is incom- 
patible with phase (Y). Therefore the phase boundary must 
be crossed and the system switches to an oscillation around 
the (Sc,Sj) focus. This then causes the system to be rein- 
jected into the (Y) phase again, producing the clockwise or- 
bit. Furthermore, regardless of initial conditions, the system 
is always reinjected at the point R thereby producing a limit 
cycle. Since Sj = Sq, the fixed points approach each other as 
P(j V /p1' V approaches 1 from above and therefore the switch- 
ing frequency increases while the amplitude decreases as con- 
firmed in (c). (b) Time series of input supply Si(t) (black) 
and catalyst Sc(t) (red) for a single process with 7 = 1. (c) 
Variation of switching frequency / against p — (Pq V /Pj' v ) 
from Eq|5| 7 = (black), 7 = 1 (red). As 7 increases the 
bifurcation point decreases and the frequency / increases. 



FigE{b)(c), but the exact bifurcation point varies. Fur- 
thermore the equivalence between switching and price 
crossing, seems to be still approximately true. How- 
ever the switching dynamics is more complex as shown 
in Fig|3^b). This is interesting since one would naively 
expect that using price information should allow the pro- 
cess to obtain more balanced input supplies and process- 
ing to occur more smoothly. However doing so moves the 
process nearer the trapping birfurcation enhancing the 
switching feedback effects. This can produce quasiperi- 
odic or weakly chaotic dynamics. 

To understand the behaviour of a multi-process econ- 
omy and the appearance of complex multi-timescale dy- 
namics, shown in Fig|21 we need to see how individual 
processes affect each other and in this respect we need 
two points, (i) Trapping bifurcation produces extremely 
different frequencies. When there are several processes 
each process has its own, now time dependent, trap- 
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FIG. 4: Hierarchical multiple timescale limit cycle attractor in 
process chain in fixed environment, (a) Diagram of a 3 process 
7 product chain. 3 products are catalysts, one in each pro- 
cess each with a fixed external supply Sf c and demand D% c , 
defining an external price p*Q = D\ c /S\c for each catalyst. 
There is one input product I with a fixed external supply Sf 
and one output product O with fixed external demand D a - 
There are two intermediate products with no external sup- 
plies and demands, (b) Time series of supplies and detail (c), 
colours are as depicted in (a). Due to price equalization Eq|5| 
in chains of processes we get p X j' y — p Q ,v = Dq/SJ for the in- 
put and output prices. Here we set external environment such 
that p x '^ >> >> V%c ~> PT V so that the 3 processes are 
all in the trapped state with very different switching frequen- 
cies, and we obtain a multiple timescale limit cycle. Since the 
3 processes are in the same chain and they are all in catalyst 
minimum state the process 3 brown catalyst oscillates around 
the output level of process 2 which is the blue catalyst which 
oscillates around the the output level of process 1 which is 
the red catalyst, to give the pleasing (but economically unde- 
sirable) effect of sewing. Here price feedback 7 = 0.9. 

ping bifurcation parameter p l {t) = p l c (t) / 'p}(t) and own 
switching frequency, determined by its external environ- 



ment, where the superscript % now refers to each process' 
'fixed point' prices given by Eq|5| The switching state 
can be turned on and off when the fixed point prices p l c (t) 
and p\(t) cross and as explained above, since the switch- 
ing frequency depends strongly on p, not on the complex 
eigenvalues of focus in the conventional way, switching 
frequencies can be extremely different. 

(ii) Unequal coupling and transfer of oscillations. We 
note that this system is an evolving system of compet- 
ing coupled oscillators that can have very different funds 
Fi{t). They are coupled only through a series of prices 
(Eqs[2EJl, so that not all processes are directly coupled. 
Furthermore the price couplings (Eq|5J| are of mean-field 
threshold type, so that each process transfers it oscilla- 
tions with different 'weights' determined by their differ- 
ent funds Fi(t) which control the size of their demands 
Dij (t) , (Eq|3J ■ Due to this the oscillations of a given pro- 
cess may dominate its input product price for example 
but hardly appear in its catalyst price, allowing multiple 
timescales and indeed trapping itself to appear. 

An example of how these (i) and (ii) combine in a 
multi-process system is the simple chain shown in Fig0] 

Summary, ft is clear that when economic processes 
with minimum condition dynamics are coupled together, 
price level crossings will induce hierarchical trapping bi- 
furcations with different switching frequencies in individ- 
ual processes. These oscillations will be transmitted to 
some prices but not others and the dynamics can easily 
become very complex as illustrated in Fig|5] which was 
produced from random networks similar to that shown 
in Fig^ This behaviour may be the fundamental origin 
of persistent cyclic non-equilibrium economic dynamics. 

The novel trapping bifurcation and switching state dy- 
namics introduced here may be relevant in other fields 
with meta-level coupling where real variable dynamics is 
synergetically mixed with boolean type logic. Economic 
agents choosing a maximum utility action, cellular re- 
action path expression interacting with genetic sorting, 
networks of molecular machines in a cell and neural de- 
cision making systems are possible examples. 
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